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I. INTRODUCTION 



The problem of finding solutions of the initial value constraints ( "initial data" ) of general 
relativity in vacuum, or coupled to matter, is of interest from both mathematical and physical 
viewpoints The mathematical questions concern existence and uniqueness of solutions for 
various classes of initial data and topology of Cauchy surfaces. The more physical questions 
concern cosmic censorship, and numerical evolution of physically relevant initial data. 

There has been much work on construction of initial data sets using analytical and 
numerical techniques, as well as combinations of both. Analytical methods initiated by 
Lichnerowitz, and further developed by York, are now standard material 0. Combinations 
of analytic and numerical techniques continue to be investigated, particularly for problems 
of physical interest. 

In vacuum, perhaps the simplest initial data is the spherically symmetric "throat" con- 
necting two asymptotically flat regions, which has a generalization to N throats with masses 
and charges [Q. There is also a wormhole solution due to Misner [^, where the wormhole is a 
handle on flat space. Other multi-black hole solutions, having time and inversion symmetry, 
have also been derived [^], with subsequent generalization by Bowen and York [^]. In more 
recent work, a variation on the theme of analytic-numerical data for black holes with 
arbitrary momenta and spin numerically is given in [0]. 

Most of the analytic solutions of the initial data problem are in vacuum. Analytic 
solutions with specific matter coupling are less studied. Such solutions are potentially useful, 
both for their intrinsic value, and for providing starting points for numerical time evolution 
schemes. As an example of the latter, analytic initial data for scalar field coupling may be 
useful for further study of the spherically symmetric collapse problem [Q, which has been a 
subject of recent interest. This collapse problem has also been studied for matter couplings 
other than the scalar field 0, but always with numerically generated initial data. 

A second reason for seeking analytic initial data sets is their potential usefulness for 
probing the cosmic censorship conjecture. One statement of this conjecture is that trapped 
surfaces in a spacetime always lie inside event horizons. This has led Penrose to suggest an 
"initial data test" for weak cosmic censorship |TD|. The test is an inequality relating the 



area of the outermost trapped surface A{S), the area of the event horizon Aeh, and the 
ADM mass M of asymptotically flat initial data, all for matter satisfying reasonable energy 
conditions. The inequality is 

A{S) <Ao< Aeh < levrM^ (1) 

where Aq is the area of a surface that encloses the outermost marginally trapped surface 
S. The first and second inequalities are reasonable from a physical viewpoint because, as 
matter collapses the region containing trapped surfaces gets larger, and in the long time 
static (or stationary) limit ^4(5*) Aeh- The third inequality allows the possibility that 
all the mass does not end up in the black hole. The content of the cosmic censorship test is 
this: If initial data exists such that the area of the outermost trapped surface is larger than 
IGvr times the square of the ADM mass of the data, then weak cosmic censorship is violated. 

Trapped surfaces may be present on an initial data surface for certain ranges of the initial 
data parameters. The boundary separating trapped and untrapped regions on the initial 
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data surface is the apparent horizon.^] On this horizon the outward expansion of hght rays 
vanishes. In terms of initial data, a closed spatial 2-surface S with normal s"' is outer trapped 
if 

{q^'-s-s'){Kab + D^s,) = 0, (2) 

where Kab is the extrinsic curvature of the initial data surface, and qab is its metric. [| This 
equation expresses the vanishing on S of the outward null expansion.^ 

Eqn. (|I]) may be viewed as providing a relation between local and global information 
in the sense that the apparent horizon is a solution of a differential equation involving the 
(local) phase space variables, whereas the ADM mass is an integral of a phase space function 
over the boundary of the initial data surface. While it is true that the apparent horizon is 
also an embedded 2-surface (and hence 'global'), its shape and size are subject to change 
due to local matter flows. 

There have been a number of successful tests of the inequality (|I]). The first were for null 
shells with flat interior, in certain special cases [0,0 . However, it has been shown recently 
that the inequality holds generally for null shells with fiat interior |13|, (see also [1^). For 



arbitrary matter satisfying the dominant energy condition, it was proven by Jang and Wald 
p!5| that time-symmetric initial data (zero extrinsic curvature) satisfies this inequality. The 
status of the inequality for the time-asymmetric cases is in general open (with the exception 
of null shells with fiat interior mentioned above ilTBlO). 



With the above motivation, it is the purpose of this paper to present classes of analytic 
initial data sets for certain matter couplings, discuss some of their properties, and use them 
as probes of Penrose's inequality. 

The matter couplings discussed are to the massless scalar field, electromagnetism, or 
Yang-Mills theory. It is based on the conformally fiat ansatz for the spatial metric 
The next section gives solutions for a class of time-symmetric cases, followed in Section 
III by time-asymmetric data for scalar field coupling. This section also contains results 
pertaining to Penrose's inequality for time-asymmetric data: the inequality is satisfied for 
all the examples considered, with the exception of cases where only the weak energy condition 
holds. 



^There is also an alternative terminology in the literature: only the outermost marginally trapped 
surface is called the apparent horizon. Here we refer to all marginally trapped surfaces as apparent 
horizons, as in Ref. |11|. 

^In terms of the gravitational momentum tt"'^ = y/q{K°'^ — K q"'^) ^ where K = K°'^qab, this equation 

is TT^^SaSb = y/qDaS"-. 

^The vanishing of the inward null expansion is obtained by changing the sign of s"" in (|2|). This 
corresponds to a white hole situation. 
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II. TIME- SYMMETRIC DATA 



Consider Einstein gravity minimally coupled to a massless scalar field with no self in- 
teraction. In the Hamiltonian formulation, the phase space variables are the canonically 
conjugate pairs {(p,P) for the scalar field, and {qab^T^""^) for the gravitational field. The 
initial value constraints are 

n = ^Gabcd^'^'r'' - ^^'^R + ^ ^qabg^^g^^ ^ (3) 

r = dbn^"" + Pq'^^dict) = 0. (4) 

where Gabcd = {gacQbd + QadQbc - QabQcd)/'^ is the DeWitt supermetric. Da is the covariant 
derivative of qab-, and ~ denotes densities of weight one. 

Consider the topology of the spatial slice S to be P? , with the rectangular coordinates 
X, y, z. The usual time-symmetric conformal ansatz for the vacuum constraints is 

= ; qab = i^Xb, (5) 

where 6ab is the flat Euclidean 3-metric. With this ansatz the diffeomorphism constraint is 
identically satsified and the scalar constraint Ti. becomes the Laplace equation. The simplest 
one throat (or Schwarzschild) solution of the Hamiltonian constraint is 

^■ = 1 + ^. (6) 

where m is a constant (the ADM mass) and = x"^ + y"^ + z"^ . To avoid the singular 
behavior at r = this solution is viewed as being on R^ — {0}. Evolution of this data gives 
the Schwarzschild solution. This result generalizes to black holes with masses mj at 
the coordinate locations on R^ — {points excised for each black hole}: 

Consider now the following ansatz for the Einstein-scalar field theory. Again with the 
spatial slice J2 ^ R^, set 

p = 0, 7^'^^ = 0, qab = ^Xb, (8) 
with arbitrary initial scalar field (j){x, y,z,t = 0). The Ricci scalar of qab is 

= -8V^i/j/ip^. (9) 

With this, the Hamiltonian constraint (|^) becomes 

+ i S'^'dacPdtcP^ ^(x, y, z) = 0. (10) 

This is the Schrodinger equation in three dimensions with potential V{x, y, z) and eigenvalue 
E related by 
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-^[E-V{x,y,z)] = -5'^' d,ci>d,4>. (11) 

Since the r.h.s. is positive or zero, the physical solutions of the Hamiltonian constraint must 
satisfy E — V > 0. This leads to an infinite class of solutions, derived from any V and 
E, if the scalar field is allowed to be arbitrary, and possibly singular on the initial data 
slice. However, since regular initial data is more physically relevant, the "wavef unctions" 
ip of interest, for example for matter collapse, will mainly be those which do not vanish 
anywhere. This criteria appears to exclude normalizable wavef unctions, where ip vanishes 
somewhere. A further restriction is that the scalar field must be real; there are E and V 
for which this is not the case for all r. These conditions are however not drawbacks to the 
analogy with the Schrodinger equation, as it is still possible to remove the normalizability 
boundary condition to get many instances of non-singular initial data sets with real 0, as 
we see below. 

(For a free massless scalar field, the Hamiltonian constraint (|10|) can also be converted 
into the Poisson equation V'^ip = — V^Vb(5"*, where the vector field Va is freely specified; the 
initial scalar field is determined from any solution ip using da4> = Vaj \f^- Also, neither this 
nor the previous case gives a linear equation for ip if the scalar field is massive or has an 
arbitrary potential term.) 

It is straightforward to obtain some simple explicit examples of regular solutions (ie. 
nowhere vanishing ?/)), and a physical picture. 

(i) Solutions in spherical symmetry: Given a radial wavefunction ilj{r) with eigenvalue E for 



a potential V{r), the initial scalar field (j){r,t = 0) from (pA]) is 



0(r, t = 0) = ^^-pfdr' ^E-V{r'). (12) 

Thus, the exact initial scalar field which solves the Hamiltonian constraint is proportional to 
phase of the WKB wavefunction for any potential V{r) and energy E, with the restriction 
that E - V{r) > to get real 0(r, t = 0). 

A particular solution in spherical symmetry: consider the scalar field pulse defined by 

r if < r < ro 

(j)[r, t = 0) = < y/2lnr if tq < r < ri (13) 
I if r > ri 

where rQ,ri are parameters giving the pulse width. One solution for ip{r) is 

{1 if < r < ro 

A/y/r + B Inr/y/r if Tq < r < n (14) 
C + D/2r ifr>ri . 

Here we have chosen the 3-metric to be flat in the "inner" region, Schwarzschild like in the 
"outer" region, and determined by the scalar field shell in the middle. The coefficients A, 
B, C and D are determined by continuity of ip{r) and its first derivative at the interfaces 
r = rQ and r = ri to be 



A = V^o{l-l Inro) ; ^ = (15) 
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Ina 



D 



'a 

tq^— Ina, 



(16) 



where a = ti/tq > gives the "width" of the pulse. The ADM mass of this solution 
M = CD is a. function of the "distance" tq of the pulse from the origin and the number a. 
It is possible to have a non-flat vacuum in the inner region by setting ijj = 1 + m/2r in the 
inner region. The data then has the additional parameter m, and is singular. 

In spherical symmetry with vanishing extrinsic curvature, the apparent horizon equation 
(0) simplifies to g°'^daRdhR = , where R = rip'^lr) is the radial coordinate of the nested 



spheres and 

gab ig 

the spacetime metric, 
this leads to the equation 



For the conformally flat ansatz for the 3-metric, 



tjj + 2rV^' = 0. 



(17) 



Thus, apparent horizons exist on the initial data surface if there are real solutions r = tah 
of this equation, with vah lying in the appropriate region(s) for each ip. For the solution 
(0), there are no apparent horizons in the region of non- vanishing scalar field because ( |T7D 
reduces to IB j ^ = 0. 

There are other examples where this is not the case, and apparent horizons are 
present. One interesting case is a solution which is asymptotically fiat without matching 
to Schwarzschild exterior. Consider 



V2C 



(18) 



where C is a constant. The corresponding solution of the Hamiltonian constraint is 

i^ = Acos (^) + B sin (^) . (19) 



This solution is asymptotically fiat. For large r, the constant A is a conformal factor which 
may be set to unity. Comparing with Schwarzshild shows that the ADM mass M of this 
data (with A = 1) is M = BC. It has a curvature singularity at = rip'^ = 0. It is 
interesting to note that there are an infinite number of zeroes of ijj and a corresponding 
number of apparent horizons near these zeroes, even though the scalar field is singular only 
at r = 0. Surprisingly, for this data there are regions of parameter space where Penrose's 
inequahty appears to be violated. One such region is at and near the point B = 0.01 and 
C = 2.00 (with A = 1), where 2M — Rah = —0.286. However, this does not provide a 
counterexample of cosmic censorship because the conformal factor ijj goes to zero outside 
the horizon |jl8|. We note that Penrose's inequality may still be tested, and is satisfied, for 



the apparent horizon that lies in the region connected to spatial infinity; ie. as one comes 
in from infinity, the horizon is encountered before any zero of ip. 

It is possible to get solutions with a finite number of apparent horizons by patching flat 
space in the inner region, from r = to some r = ri. Set ri = 1, and = and ip = 1 for 
< r < 1. Then, with (p as in (|T^) for r > 1, the solution is 



V^(r) = A{C) 



cos 



2r 



tan — sm — 
2/ V2r 



for r > 1, 



(20) 
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where 



A{C) = [cos(C/2) + tan(C/2) sm{C/2)]-\ (21) 
The ADM mass of this one parameter (C) configuration is 

M = A^iC) C tan (^^j . (22) 

Now, depending on the value of C, there are zero or a finite number of apparent horizons in 
the matter region. It is also possible to construct a solution with a finite number of horizons 
by patching Schwarzschild rather than fiat data in the inner region. It is again possible 
to choose parameters such that a horizon is encountered before a zero of tp, as r is varied 
inward from infinity. 

(ii) Scalar field wormhole: Away from spherical symmetry, the situation is virtually un- 
changed. Consider perturbing Misner's wormhole solution |Q by the presence of a compact 
scalar field pulse. The mathematical problem is one of matching an 'interior' vacuum solu- 
tion to a pulse solution, which in turn is matched to an 'exterior' vacuum solution. (Vacuum 
solution means a solution of the constraint equations without matter sources.) The former 
and latter are just Misner's ip, with different parameters in each region to account for the 
mass in the scalar field. The solution in the intermediate region depends on the scalar field 
pulse, which determines the potential in the Schrodinger equation analogy. The details are 
a matching problem not unlike the one in spherical symmetry considered above. Evolu- 
tion of such initial data would be of interest for seeing how presence of matter affects the 



gravitational radiation from the two black hole problem in the so called "close" limit [|T9 . 

Electromagnetism: There are similar results for coupling to electromagnetism, with phase 
space variables (A„,E") and Hamiltonian density {E'^E'' + B''B^)qab/2^, (5° = e'''"'dbAc). 
Using the same conformal ansatz for the metric and conjugate momentum, qab = ip'^Sab and 
tt"^ = 0, set 



= e" = (23) 

where e°'{x,y,z) and b°'{x,y,z) are arbitrarily specified. Then the Hamiltonian constraint- 
becomes 



'^' + ^{e^eXb + b''b%b) 
lb 



^ = 0, (24) 



identical in form to the scalar field coupling case (0). Thus, given the freely specified fields 
e"' and 6", one solves the Hamiltonian constraint for ip, and then determines the physical 
electric and magnetic fields via (|23|). Note that this ansatz is not the same as the one 



given by Misner and Wheeler |0, which gives charged black hole data; their metric ansatz 
is Qab = ("0^ — x'^)'^^ab, from which the scalar constraint gives Laplace equations for ip and x- 
As for the scalar field case, there is an alternative specification of the electric and magnetic 
fields which converts the Hamiltonian constraint into the Poisson equation with source given 
implicitly by these fields. This arises by setting 

E" = ^^^"^e" B" = ij^/^b", (25) 
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and gives 



= -T^(^"^'^-b + b%X,) (26) 
lb 

for the Hamiltonian constraint. 

Of course, neither (^) or ( pSj ) solve the diffeomorphism constraint unless either the 
electric field or vector potential Aa is set to zero. Furthermore, there is also the Gauss law 
for the electric field. One obvious solution to all the initial value constraints is obtained by 
setting E°- = 0, and solving the scalar constraint or (|26|) for purely magnetic initial 
matter. In this context, one can again consider special cases, such as spherical or toroidal 
symmetry, or a wormhole. 

Purely electric solutions, other than the monopole, are more difficult to find without 
imposing further symmetries, since it is difficult to simultaneously solve the vacuum Gauss 
law constraint anc? obtain a linear equation for the scalar constraint. However, this is possible 
if symmetries are imposed. As an example consider the case of one translational symmetry 



and seek solutions ip = ijj{x,y). For two arbitrary functions u{x,y) and v{x,y), set ^0 

~a ^ ^abcQ^^g^^^ (27) 

where e"'"^ is the metric independent Levi-Civita tensor density. Then the electric field given 
by ( P^ satisfies Gauss's law 

daE^ = daiil^'^e'^^^dbudcv) = 2ilje^^''da^dbudcV = (28) 

because u, v and ip depend only on two coordinates. The electric field lines are tangent to 
the curves defined by the intersection of the surfaces u = constant and v = constant. In the 
present case, the only non- vanishing component of E"" is E^{x,y). 

Finally, we point out a class of solutions for which both initial electric and magnetic fields 
are non-zero. This is again in the context of one translational symmetry. Given arbitrary 
u{x,y), v{x,y) and b"'{x,y), define e°'{x,y) according to (P?]). As before, since the solutions 
of the Hamiltonian constraint are ip = ilj{x,y), Gauss's law is automatically satisfied. The 
diffeomorphism constraint 

dar' + E\daA, - dtAa) = (29) 

reduces to 

iP^lahc&B' = 0. (30) 

Now, since E"- has only a component, restricting B"^ to have only a z— component solves 
this constraint. Furthermore, daB"^ = dz{ip'^{x,y)b^{x,y)) = so there are no magnetic 
sources. 

This procedure gives solutions of all the initial value constraints with electromagnetic 
coupling; the data is characterized by two arbitrary functions of two coordinates. As there 
are no electromagnetic sources, this class of solutions may be called "geon" data. This result 
has extensions to any 3-space with one Killing symmetry. 
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Yang-Mills theory: The phase space variables are (A^, where i (= 1 ■ ■ ■ A^^ — 1) is 
the SU(N)) Lie algebra index. The Yang-Mills contribution to the Hamiltonian constraint 
is 

H = ^{E-'E^ + B'''B'^)qabh, , (31) 

where kij is the Cartan metric, and 5'^* = e'^^^F^. (Fj^ = daAl - d^Ai + C'ji^A{A^ and C'j,^ 
are the structure constants of the gauge group). The ansatz that reduces the Hamiltonian 
constraint to the Schrodinger equation is similar to (^), namely 

Thus, given arbitrary e*^* and 6"*, finding the conformal factor ip is again an elementary 
problem. However, we must also find solutions of the diffeomorphism constraint and the 
non-abelian Gauss law 

DaE"' = daE^' + C'j^A'^E''^ = 0. (33) 

One class of solutions of all the constraints is obtained by starting with = 0. This 
immediately solves the diffeomorphism constraint (because vr"^ = 0), and the Gauss law 
becomes dai'ip'^e"'^) = 0. Now set 

= C'ji^e'''"'dbU^d^v'' (34) 

where and are 2{N'^ — 1) arbitrary functions. This form does give a "full" solution of 
the Gauss law for this case (ie. with A^ = 0). This is because its solutions, without any 
redundancy, are parametrized by 2{N'^ — 1) functions for SU{N), since there are A^^ — 1 
constraints for the 3(A^^ — 1) E°'^ fields in three spatial dimensions. With this ansatz Gauss 
law becomes 

9„(7/>2g«) = 2^C'ji^e''^^datpdbU^dcv\ (35) 

Then, as for the abelian case, there are solutions to all the constraints if the problem is 
reduced so that all variables depend on two coordinates. 

It is more difficult to find solutions for which both electric and magnetic fields are non- 
vanishing. This is due to the term quadratic in the phase space variables in the Gauss 
law. To see this, suppose we give an ansatz specified by an arbitrary gauge field a\ whose 
magnetic field is 6*^% and with e"* given by 

g'^^(x) = J2 e^^'dbU^'d^v'' Tr[f/[a„, 7(m^, ?;^)]t*](x) , (36) 

a 

where and v'^ are a set of scalars (labelled by indices a, b,), r* are matrix generators of 
the group, and 

U[Aa,j{u,v)]{^) = Pexp / dsAi{j{s)y (37) 

J-y 

is the holonomy of the gauge field along the loop 'y{u'^,v'^), with base point x, determined 
by the intersection of surfaces u'^ = constant and v'^ = constant. This is a non-abelian 
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generalization of (pT]) pO[ . It is straightforward to check that it satisfies the Gauss law for 
the fields (e'^*,a^). However, because of the conformal factor ip in (|32|) , it does not satisfy 
the Gauss law ( p3D for the physical fields {E°-\ A^). Indeed, from a solution ip of the scalar 
constraint obtained via (021), we find for (|36|) that 



DaE'^' = Daii,^^') = ij'^f^Mi - + 2ije^'dai^ ^ 0, (38) 

where is a gauge field associated with the physical magnetic field 5"* = ip'^W'^. 

One way to solve (^) is to assume, at the outset, that and are parallel in the 
internal direction, ie. that e"* = e"t* and = a^t* for a fixed Lie algebra vector f. This 
means that [a^, Ob] = 0, and and ^a('?/', Aq) are also parallel in the internal direction. With 
these conditions on the initial configuration, the first term in ([38| ) vanishes. The second term 
vanishes if we assume, as before, that all fields depend on only two of the spatial coordinates 
{x,y), ie. that there is a translation Killing symmetry. This means that the YM electric 
field from (^) has only a 2;— component. 

The diffeomorphism constraint 

eabcE'^B'^^h, = (39) 

is satisfied if the YM magnetic field is also restricted to have only a z— component. Finally, 
note that DaB""^ = d^B'''' + [A^yB^f = because -B""* is a function only of x,y, and the 
commutator vanishes by the ansatz used to solve the Gauss law. Thus there are no magnetic 
sources. 

In this way one can obtain solutions to all the constraints for the time symmetric Yang- 
Mills problem, with the conformal ansatz for the metric. Such solutions are parametrized 
by functions of two variables and v^. Because of the way the Gauss law is solved, the 
Yang- Mills matter for this data are arranged into "non-abelian lines of force." This data is 
therefore a direct generalization of the electromagnetic data given above. 

Charged scalar field: For gravitational coupling to both electromagnetism and a charged 
scalar field, the Hamiltonian constraint has the additional term 

^q''''Va(j)V,(j)* + PP*/2^, (40) 

where T>a = da+Aa and * denotes complex conjugation. With the time-symmetric conformal 
ansatz, and P = = Aa, the Hamiltonian constraint is 

8VV + + e^e'6ab) ^ = 0, (41) 

which must be solved with the Gauss law daE'^ = Pcf)* + P*(j). 

A large class of solutions of this system can be obtained when there is one translational 
symmetry. As before set e*^ = e°'^'^dbudcV, where u and v are arbitrary functions of the 
coordinates x, y, z. (p is also arbitrarily. Then the linear equation (^) can be solved for 
With a solution ijj, the Gauss law for E"- = ip'^e"' (with P = 0) becomes 

e'^da^ = i'^'^dbUdcVda^ = 0. (42) 

Thus, if all the functions u, v, ip depend on only two coordinates, the Gauss law is identically 
solved. 
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III. TIME-ASYMMETRIC DATA 



So far we have restricted attention to time-symmetric (ie. tt"'' = 0) initial conditions. 
We now consider some cases which relax this condition, but still restrict attention to the 
conformally flat form of the three-metric. The goal of the ansatze is to obtain more general 
initial data by attempting to convert the Hamiltonian constraint into a solvable equation. 
This is useful not only for obtaining more general classes of solutions, but also for addressing 
the cosmic censorship conjecture via the initial data test (|1|); as noted in the introduction, 
the status of this initial data test is open for the general time-asymmetric case. 

If Tc""^ 7^ 0, the vector (or spatial diffeomorphism) constraint must also be solved. It 
is obvious (and well known) that in vacuum an easy way to solve this constraint is to set 
^ab _ j\qab ^ where A is a constant. When matter is present, it is possible to obtain initial 
data by generalizations of this, as we now describe. 

Consider the scalar field case with the spherically symmetric ansatz 

qab = ij\r)5ab n'^' = iJ-'P{r) [an'^n' + (36'^') , (43) 

where P{r) is the scalar field momentum density, = x"/r is the unit radial vector, and a 
and /3 are constants. The diffeomorphism and Hamiltonian constraint become 

{a + P)dr{'ip-'P) + + 2p-^Pdr(P = 0, (44) 



— - a/3 - ^ + 1 ) ij-^P'' + + i?da^ 



0. 



(45) 



The diffeomorphism constraint gives 



P 



r 



exp 



a + (3 



(46) 



One class of solutions is obtained by arranging cancellation of the gravitational and scalar 
momenta. Thus, with P arbitrary, set the coefficient of P^ to zero. This gives 



a 



/3 l±2Wl- 



1 



2/?2 



(47) 



Now, if the scalar field is chosen to fall off sufficiently fast at large r, qab falls off like 
Schwarzshild at spatial infinity. However, from (^) it is evident that the leading order term 
of tt"^ cannot be 0(l/r^); the falloff is slower than this. For the special case = 1 (flat 



slice), the required falloff is r 



-3/2 



^see below). Even this violates (^T]). This means that 



this data is not asymptotically fiat. 

There is a generalization of this data obtained by requiring that d, 
P, again with ip = 1 in (^). This requires (p ' 
this also does not give asymtotically fiat data. 



be proportional to 
Inr, which means that P ~ r^^. Therefore 
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In order to find asymptotically fiat data, let us consider the general spherically symmetric 
ansatz with fiat spatial metric ie. ip = 1, or qah = ^ab- The general form of vf"'' is 

= /(r)nV + g{r)5^\ (48) 

The diffeomorphism constraint, with matter current J*^ = 0, is 

/(r) = -^,Jdr {r%.g) + ^, (49) 

(c = constant) and the Hamiltonian constraint is 

^"'^ab - \{K? + P = \U + 9)U~^9)+P = ^, (50) 

where p is the matter energy density. 

In vacuum (p = 0), the solution of the constraint equations may be found by taking 
g( = r" in the above equations. The solution is 

^'^'=^(372^ + 5'^''), (51) 

where C is a constant. Since this is vacuum data is spherical symmetry, it must be data 
for the Schwarzschild metric. Its mass M may be determined in terms of C by locating the 
horizon on the initial data surface. With = in the apparent horizon eqn. (^, and the 
horizon radius fixed to be 2M, eqn. gives 

C = (52) 
The corresponding space time metric is [BlI] 



The fiat slice data for Schwarzschild is unusual in that the normal to the spatial slice is 
not perpendicular to the normal of the timelike boundary. Because of this, the mass formula 



is different from the standard ADM integral in that it involves the extrinsic curvature ^2 
The above approach provides a direct way of obtaining fiat slice data for the Schwarzschild 
metric. 

From (^) it is evident that even with matter, fiat slice data must have tt"'' falling off 



at least as fast as r to maintain asymptotic flatness. Furthermore the positive energy 
condition p > requires from (|50D that 

if + 9){f-^9)<0 (54) 

for all r. It is difficult to flnd analytical data explicitly that satisfies both these conditions, 
and (ff^), without patching to Schwarzschild exterior. Consider for example the electric 



"^I thank Ted Jacobson for suggesting the use of flat slices. 
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monopole. The solution of daE"" = is i?'^ = Qn°- jr^ (with the point r = removed), and 
the magnetic field F^h = 0. Then the condition on /(r) and g{r) from the diffeomorphism 
constraint is (H9|), and that from the Hamiltonian constraint is 



Q 



2 



f - 2fg -3g' + ^ = 0. (55) 
These two conditions give a single equation for /(r) (or g{r)): 

f' + -+'^X' = 0, (56) 
r 2 

where e = ±1 and X = ^/Af^~+'3E^, {E = \E°-\). The same equation applies to scalar field 
coupling or other matter coupling with the appropriate replacement of energy density. 

It is possible to find solutions of this equation numerically. The large r behaviour is 
~ i'-3/2 expected, if the matter falls faster than r~^/^ at large r. For e = +1, / and g are 
positive for large r, so the mass of the data is positive, and there is an apparent horizon. 
For e = — 1 on the other hand, /, g, and hence the mass are all negative, and there is 
no apparent horizon. This is surprising because the matter satsifies the dominant energy 
condition. It suggests a naked singularity like the negative mass Schwarzschild case. 

It is easier to find analytic solutions, with matter satisfying p > 0, if the data is taken 
to have Schwarzschild exterior. An example is provided by taking the interior and exterior 
jfab gjygj^ g = (y -\- p-f, (which implies / = — /3r/3 (1^)), and (|5l|) respectively. Matching 
jj-ab ^ = I gives 





(57) 

and interior energy density 

p = 15M(1 -r)(5-3r), r < 1. (58) 

This energy density may be interpreted as arising from a scalar field by setting p = and 
= 0, or from electromagnetism by setting p = E°-E^6ab, B'^ = 0, with the Gauss law solved 
as in the last section. 

The apparent horizon equation for ( ^HD reduces to /(r) +g{r) = 2/r, and for this example 
gives 




This shows that an apparent horizon forms for M = 2(12/25)^ with non-zero radius r = 5/6. 
Using (p9[), it is possible to show that Penrose's inequality is satisfied for all M. Other 
examples, with / and g polynomials also satisfy the inequality. 

If the condition ( P^D relating /(r) and g{r) is not imposed, then J" ^ 0; the only 
restrictions on these functions are now those arising from energy conditions. An example of 
fiat slice asymptotically fiat data with J" 7^ is obtained by setting 
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m = ^, 9ir) = (60) 

Then the diffeomorphism and Hamiltonian constraints give 

M_ (a. + 3,a)(,a-a) 

respectively. This data satisfies the weak energy condition for all r if a and f3 are such that 
p > 0, but violates the dominant energy condition p > {J^Jaf^'^ for r > (a + 3/3)(/5-a)/4/3. 
The mass of the data is zero because the large r falloff of tt"* is faster than the r required 
to get a non-zero mass from the relevant surface integral |]22|. (Note that it is possible to 



get a zero mass even if there is a non-zero energy density, provided the energy density gives 
a metric whose leading order behavior is qab ~ ^ab + 0(l/r^); ie. there is no 1/r term to 
be captured by the surface integral. Conversely, a form of qab and vr"'' such that the ADM 
mass is manifestly zero, may be used to deduce an energy density and matter current. This 
is exactly what is done above for fiat slice data. Although this seems counterintuitive, it 
is possible: a metric with zero ADM mass but non-zero energy density is the Reissner- 
Nordstrom metric with M = and electric charge Q 7^ 0; while there are no horizons in this 
case, there are for the above fiat slice example.) 

The apparent horizon equation is a — (3 = 2r, showing that there are no horizons for 
f] > a > 0. Therefore Penrose's inequality holds even though the dominant energy condition 
does not. However there are other examples where this is not the case. Consider e.g. 
f3 = —a < 0, for which p = 2a^/r^; there is a horizon at r = a. Thus, the inequality can be 
violated if only the weak energy condition holds. 

More generally, from ( |48D and (0), the energy condition p > {J"'JaY^'^ is 

\{f + 9){'i9 -f)=c(^f + ^ + g'^> 0, (62) 

for constant c > 1. Simple ansatze such as f ^ g, or f = gh for some function h, lead 
for several analytically solvable cases, to data which is not asymptotically fiat. However, 
for any given asymtotically fiat form for / or (7, it is relatively straightforward to study the 
resulting ordinary differential equation numerically to probe Penrose's inequality. 



IV. SUMMARY 

We have given a number of ansatze for solving the initial value constraints of general 
relativity with matter couplings. These include the massless scalar field, electromagnetic 
and Yang-Mills fields. Large classes of solutions are obtained in each case, both for time- 
symmetric and asymmetric situations. For the scalar field case, explicit initial data sets 
are given for pulses. These may be used as staring points for numerical integration. For 
electromagnetic and Yang-Mills fields, the data may be interpreted as the initial geometry 
due to fiux fines of the electric field, because of the way the Gauss law is solved in each of 
these cases. These results go beyond the multi-point sources considered in earlier works. 

For the time-asymmetric cases considered, Penrose's inequality holds if matter satisfies 
the dominant energy condition. However, violations of the inequality can occur if only the 
weak energy condition holds. 
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The solutions given in this paper may be used as a possible starting point for extensions 
away from spherical symmetry. This would be useful not just for the data that can be 
obtained, but also for providing new tests of cosmic censorship via Penrose's inequality. 
This is a potentially important direction because almost all results pertaining to cosmic 
censorship are in spherical symmetry. 

Several generalizations of the ansatze given here arise if a fixed direction S"" is specified. 
Then n"-'^ can be constructed out of 5", n"', and the metric. One may even take S"" to be a 
non-constant divergence free vector field specified similarly to the electric field which solves 
Gauss's law. These and similar extensions are presently being studied. 

Acknowledgement: I would like to thank Ted Jacobson and Bill Unruh for discussions. This 
work was supported by the Natural Science and Engineering Research Council of Canada. 

Note added: After this work was submitted and posted, I learned that Penrose's inequality 
has been previously studied in spherical symmetry: it was proved under certain conditions 
in [^, and more generally in 1^^. I thank E. Malec and S. Hayward for bringing these 
references to my attention. 
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